Extended and localized phonons, free electrons, and diffusive states in 

disordered lattice models 
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In this paper we propose that phonons, free diffusive electrons, and diffusion in two-and three-dimensional disordered lattice 
problems belong to a different class of localization than bonded electrons. This is manifested by three effects that can be 
observed numerically. First, there are extended states even at two dimensions, whereas there are no extended states in the 
usual electronic models. Second, the correlation length does not diverge at the mobility edges in three dimensions, and finally, 
the participation ratio of the extended states, decays to zero at this edge. This indicates zero electronic conductivity, in the 
extended region near the mobility edge. We show that low energy modes for these models can either have diverging localization 
lengths or are extended. 

PACS numbers: 63.20 localized modes, 63.50 Vibrational States in Disordered systems, 73.23 Theories and Modes 
for Localized States 



It is well known that disorder in the properties of physi- 
cal systems can induce localization effects on the relevant 
variables. A classic example for such phenomena is an 
electron in a disordered potential [1-10]. Specifically, an 
electron that obeys the Schrodinger equation, with spa- 
tial disorder in the potential, is known to be localized for 
a strong enough disorder, which affects the conductivity 
of an electronic system and its dependence on the system 
size. This subject has received enormous attention dur- 
ing the last three decades and is described in a numerous 
books and reviews [1-10]. 

In one dimension and for any kind of disorder, all elec- 
tronic states are localized, i.e. the wave-functions decay 
at a large distance as exp(— |r — ro|/£), where £ is the lo- 
calization length, and r and r$ are local positions. This 
leads also to exponential localization of the electronic 
conductivity. In contrast, in two dimensions the situa- 
tion is marginal. For any disorder the wave functions 
form a multi-fractal set in the system. In three dimen- 
sions there are localized and extended states separated by 
a mobility edge at an energy E c . The correlation length 
diverges in the localized region as £ ~ (E — E c )~ u , and 
the conductivity is proportional to (E — E c ) v in the ex- 
tended range ||||. Classical waves (photons and waves), 
are thought to behave similarly (see in p|]lO|]). 

Specifically, there are three equivalent physical systems 
in which the disorder appears in a random Laplacian form 
and can cause dramatic effects. The first is vibration 
problems (or phonons) in disordered media. Examples 
for such problems are disordered solids, glasses, gran- 
ular materials, packings, colloids, and polymers, where 
the geometrical disorder, connectivity and fluctuations 
in the elastic constants play a crucial role in determining 
of the dynamics. The second is free electrons on a lat- 
tice, with randomness in their real, positive transmission 
coefficients. The third is the random diffusion problem. 



The dynamical structure of all these problems is simi- 
lar though the time-dependence, statistics, and physical 
properties are very different. 

In the one-dimensional case, it is well known that 
phonons are localized [ |T2||r^ , |T4| ] . However, the local- 
ization length of these systems is scaled, as the inverse 
square of the frequency. So at low frequencies, we observe 
phonon-like modes with a localization length that is much 
longer than their wave length, and with a low-frequency 
distribution function similar to the ordered one. In a 
few theoretical papers published on the subject |p^ , |l5| , 
it was proposed that the situation in higher dimensions 
is analogous to the BE situation. We claim in this paper 
that there is sufficient numerical evidence that this is 
not the case. 

We begin by demonstrating that the dynamical ma- 
trices of the three cases are equivalent within certain 
parameter limits. Then we analyze numerically a lat- 
tice problem that is similar to the models studied in the 
electronic case [1-10]. We used increasing lattice sizes 
to get information about the nature of the eigenstates 
as the system size is increased. The participation ratio 
(PR) of the states was used to measure the degree of 
localization. A localization transition between extended 
and localized states was observed in the spectra in two 
and three dimensions. However, the three-dimensional 
transitions are different from the bonded electron (BE) 
transition[l-10]. There is no divergence of the correlation 
length near the mobility edge E c , (i.e. v = 0), and the 
volume of extended states (defined by the participation 
ratio) converges to zero at the mobility edge. This im- 
plies zero electronic conductivity near the mobility edge. 
Moreover, this can also have a major effect also on the 
heat conductivity for phonons in glasses where there are 
well-known anomalies (a plateau in heat conductivity in 
medium temperatures, see for example in WW). Further- 
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more, in two dimensions there are three kinds of modes: 
extended, localized, and multi-fractal modes separated 
by mobility edges. We propose that the difference be- 
tween these problems and the BEs is because that po- 
tential terms destroy the low-energy free structure of the 
free electrons and phonons. 

Let us define the equivalent models that we study. The 
first is the discrete random elastic Hamiltonian 

H=\^K ij {u i -u j ) 2 , (1) 

ij 

where Ui represents scalar elastic movements and K+j rep- 
resents the elastic constants. Though the us are generally 
vectors, we will limit our discussion here to scalars. There 
are two obvious constraints on the parameters. The first 
is the symmetry Kij = Kji- The second is that for all re- 
alizations of Ui, the Hamiltonian is non- negative. This is 
a non-local requirement for the realizations of the elastic 
constants that can be trivially satisfied by non-negative 
which is used. The single atom dynamics is defined 

by 

miUi = y^Kjj(ui -Uj), (2) 

ij 

where m, is the particle mass. Assuming that it, = 
exp(i(jjt)Ui{uj) we obtain an eigenvalue problem 

-m,iUj 2 Ui = } Kjj(Ui - Uj). (3) 

ij 

Free random diffusion has a similar dynamical form. 
We define densities as pi, and random diffusivities Dij 
between nearby sites. The currents will be Jij = Dij(pi — 
Pj) and we obtain a random diffusion equation 

Pi D ij(Pi ~Pj)' ( 4 ) 

ij 

Defining pi(t, p) — exp(—pt)pi(p) we obtain a similar 
eigenvalue problem with p = uj 2 . Since negative densities 
in p are not allowed, > 0. 

A free randomly moving quantum particle is described 
by 

!i=^M 13 (^-^) (5) 

ij 

for arbitrary hermitian My. For real and non- negative 
Af's we again obtain a similar equation, for ip = 
cxp(iEt)-ipi(E). These linear models have the same lin- 
ear dynamical structure for the phonon case when the 
masses are the same, and in the diffusive case where the 
diffusion constants are symmetric. However, the time 
dependence, the statistics and other properties are dif- 
ferent. In this paper we use, unless stated otherwise, the 
phonon terminology 



Let us concentrate on an equal mass (rrii = 1) model in 
a cubic lattice, in d dimensions, with a unit distance and 
size L. We consider the nearest neighbors' interactions, 
with two elastic constants K\ and Ki, with a ratio a — 
K% / K\ , between them and a probability p to find K\ and 
I — p for K%. Since a < I is the same as a > I, we will 
only discuss this limit. The normalization K\ = I defines 
a scale of frequencies. This model was suggested in 
and was simulated in |l7j for small as. 

We can calculate the density of states N(uj) and the 
eigenvectors u r {uj). The participation ratio is our main 
measure for the nature of the states. It is defined as 

h(u,) = C£<(u>))- 1 , (6) 

r 

where the summation encompasses all lattice sites. If a 
state is localized, the non-normalized PR is not depen- 
dent on the system size. However for extended eigen- 
states the PR is scaled as L d || . If the states are multi- 
fractal, as is reported for electronic states in two dimen- 
sions, the PR is scaled as with (3 < d. Another mea- 
sure of the localization is the behavior of the correlation 
length and the conductivity near the mobility edge. 

There are two simple limits for this model. Namely, 
a = 1, which is the ordered lattice limit, and a = 0, 
which is the percolation case. With a = 1, all states 
are extended and ordered. The percolation limit how- 
ever, is more involved. Above the percolation threshold 
it is expected to see phonons on a scale much larger than 
the percolation correlation length |p^| , and fractal exci- 
tations called fractons, below this length |f6| . Below the 
percolation threshold, the system would be completely 
localized. The two-and three-dimensional percolation 
cases were recently studied in [Q. Different behavior is 
observed for a > 0. Fractons with phonons were observed 
for the case of extremely small a |f7| . 

The eigenvalue problem is treated numerically using 
the LAPACK package with real symmetric banded sys- 
tems for the equal mass model. For the variable mass 
model we used a special routine for real-banded systems. 
Specifically, we used a cubic system of size L, where the 
distance between particles is unity. Most of our simula- 
tions were done for free boundary conditions. However, 
we verified that the periodic boundary conditions did not 
change them. We computed a set of eigenvalues LOi, and 
eigenvectors denoted as u r (u)i), where r is the lattice po- 
sition. 

We simulated the one dimensional case to check our 
method. The numerical results here are completely in 
accord with the theoretical predictions. We found that 
the PR scales in the low frequency regime as w~ 2 , with 
a lower frequency cutoff whose size is reduced with an 
increase in system size. 
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is shown in Fig. |[ where the scaling exponent [3 of I4 
versus L is presented as a function of to. Five scaling 
regimes of (3 are visible in the curve. There is a marked 
transition between the extended state (3 = 2 and the 
state with /3 = 1, within our numerical accuracy. Thus, 
there is a sharp mobility gap between the multi-fractal 
and extended states in this model. 

The the same behavior was observed for a very wide set 
of parameters for a < 0.5, and the probabilities between 
0.2 and 0.8. A phase diagram of those limits will be 
published in a forthcoming paper. The same effects were 
also observed for continuous distributions of the elastic 
constants. 



FIG. 1. The density of states JV(oj) for p = 0.4 and 
a = 0.1. Note the two distinct peaks that are typical to 
the two-component system. 
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In the two-dimensional case, we observed that for 
small values of a, localization effects are observed over 
a wide range of ps. Figures 1-3 present numerical results 
from a system with parameters p = 0.4 and a = 0.1. 
Figure || shows the density of states. 
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FIG. 2. The participation ratio I4 in a two-dimensional 
system for p = 0.4, a = 0.1 and various system sizes. The 
curves are averaged over equal intervals of u) of size 0.05 to re- 
duce the I4 fluctuations. The extended modes are seen below 
uj ~ 0.6. 

In a two-component system we observe two distinct 
peaks in the distribution function |]l5| . At low frequen- 
cies the density of states is linear in so the states are 
phonon-like. Figure ^ presents I±(uS) as a function of 
to. At low frequencies there are extended modes whose 
PR scales with the system size. There is a well-defined 
mobility gap at w c , separating the extended and multi- 
fractal states. The PR below it is l£ ~ (u c — lo)L 2 . To 
analyze these states we calculated the exponent (3. This 



CO 



FIG. 3. The scaling exponent /3 as a function of uj. The 
curve is averaged over intervals of frequencies to to reduce the 
noise in I4. 

The same kind of model was also simulated in three 
dimensions, where we observed a localization transition 
(see Fig. ||). The scaling of the extended eigenstates is 
again if ~ \to — lo c \L 3 . Above the mobility edge all states 
are localized. There is no divergence of the correlation 
length near the transition £ ~ (to — lo c )° (y — 0). There- 
fore this localization transition is different from that of 
BE. One would expect to see no power law contribution 
to the zero conductivity of free electrons near the E c . 
This result is related to the decay of the PR near E c . 
States with a low PR cannot transfer current effectively, 
because according to the Kubo-Greenwood formula, the 
conductivity is related to overlaps between wave func- 
tions. This is also a clear indication for zero conductivity 
in the two-dimensional case (see ||). Indeed, calcula- 
tions of conductivity in such systems will reveal a strong 
dependence of the conductivity on E but will probably 
show only slight or no dependence of the conductivity on 
the system size. Moreover, this will also have an impor- 
tant effects on the heat conductivity in glasses [0 . 

A common feature of all these models is the existence a 
minimal zero energy state. This is a sign for the existence 
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of 'free' states in a low-frequency region. To show that 
the localization length should diverge, let us consider the 
random quenched model of diffusion with minimal diffu- 
sivity. Any local density will diffuse to infinity with a 
minimal effective diffusion constant (above the minimal 
value of D). If all states are localized, no diffusion to 
infinity would be possible, since the initial condition is 
expanded by eigenstates of the problem. Therefore the 
correlation length should either diverge at low frequen- 
cies or the system will be extended. 
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FIG. 4. The normalized PR in three dimensions for p = 0.8, 
a = 0.1. The solid line is for L = 8, dotted for L = 16, and 
dashed L = 22. The curves are averaged in frequency to 
reduce the fluctuations. An average of 6 runs was used for 
L — 8. The extended modes are observed below uj ~ 1. 

Let us now make a detailed comparison using the the- 
ory. Note that both theoretical works Jl4|,[lJI are based 
on approximations which also fail in the electronic case 
||. One can compare a numerical simulation with the 
theory even for a system with a limited scale. The basic 
concept used in our simulation is the following: if there 
is an increase in the system size, localization lengths that 
do not depend on the system size will appear in the simu- 
lation when the system scale becomes large enough. This 
is what is observed in numerics in one dimension. As the 
system size is increased, there is a decreased lower zone 
of frequency where PRs depend on system size, and an 
increased range where the PRs fit the theoretical pre- 
dictions. Interestingly, this kind of effect does not oc- 
cur in two and three dimensions. Specifically, there is a 
clear localization transition in two dimensions, but the 
states above it, even the ones with the smallest normal- 
ized PRs continue to scale with the system size. A two- 
dimensional scaling length £ ~ exp(— ui~ 2 ) was predicted 
in juj. This indicates that above a cutoff given by L ~ £, 
all states would scale with the system size and that there 
would be a shift in the cutoff when L is increased. Both 
effects were not observed in our numerics. In three 



dimensions it was claimed that a blowup of the correla- 
tion exponent occurred with a similar exponent in the 
electronic case [Q. This was also not observed. The 
results of |l^] for three dimensions can also be tested, 
since they are detailed, and there is no agreement. Note 
that for larger system sizes, there might be difficulties in 
scalings; this is the common danger of numerical work. 

All the effects discussed in this paper seem to be con- 
served for different masses and also for different continu- 
ous probability distributions for the Ks. 

The findings presented here raised the following ques- 
tions that are currently under study. What is the theory 
of the observed transition? What is the effect of adding 
the mass and vector properties to the phonon transla- 
tions? What is the nature of the transition between 
bonded and free electrons? How does the electrical and 
heat conductivity depend on the energy? What is the 
connection of those effects to real glasses? 

We thank T. Kustanovich, R. Zeitak, O. Gat for useful 
discussions of the subject. 
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